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ABSTRACT 



This paper presents a method for Acoustic Emission (AE) source localization in isotropic plate-like structures based on 
the Extended Kalman Filter (EKF). The accuracy of the traditional triangulation methods depends on the time of flight 
(TOF) measurements and on the group velocity assumption so uncertainties in both should be taken into account and 
filtered out. An algorithm based on the Extended Kalman Filter (EKF), capable of filtering out these uncertainties, has 
been developed for the estimation: 1) the AE source location and 2) the wave velocity. Experimental tests have been 
carried out on an aluminum plate to show accuracy and robustness of the proposed approach. 

Keywords: Acoustic Emission (AE), Triangulation, Time of Flight (TOF), Lamb wave, Extended Kalman Filter (EKF) 



1. INTRODUCTION 

In the past two decades, significant efforts have been made towards the development of integrated structural health 
monitoring (SHM) systems in order to increase safety of aircraft systems and reduce life cycle costs 1 . Damage location 
using Lamb waves in plate like structures, can be achieved under either an "active" or a "passive" mode. Under the 
"active" mode, diffractions of piezoelectrically-actuated waves can be used to locate existing damage in a post-mortem 
mode 2 ' 3 . Under the "passive" mode, growing damage (e.g. fatigue cracks), or sudden impacts can be located by 
monitoring acoustic emissions in a real-time mode 4,5 ' 6 ' 7,8 . This paper focuses on the "passive" monitoring technique 
using a sparse array of piezoelectric transducers to locate the point of AE source. Traditionally, damage or impact 
location is based on time-of-flight (TOF) triangulation of wave measurements taken at multiple receiving points. This 
technique works very well when the wave velocity (V g ) in the test material and the arrival time (?,) of the signal at all 
three sensor locations are known. The damage or impact location is identified by drawing three circles of radii (R,), 
whose centers coincide with the three sensor locations. The radius Ri is obtained by multiplying the time of arrival of 
the signal (?,) with the wave velocity (V g ). The intersection point of these three circles is the damage location (figure 1). 
However, the exact time of arrival of the signal may not be known either because of the distortion of the signal due to 
dispersion phenomena or because the instrument noise hiding the exact time of arrival of the signal 4 . In addition, 
environmental interference, such as temperature changes, significantly affects the TOF measurements of acoustic waves 
and wave velocity 9 . As a result, rather than the damage being located at a single point at the intersection of the circles, 
the damage can be located anywhere in the dark overlapped region at the intersection of the rings around the sensors 
locations (figure 1). The ring's width represents the uncertainty in the measured distance as a result of TOF and/or wave 
velocity uncertainties. In general, the lower the signal-to-noise ratio (SNR) the higher is the ring's width (i.e. higher 
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uncertainty in TOF measurements). As a result of mentioned uncertainties the AE location, TOF and wave velocity are 
random variables and the procedure of AE source localization can be treated as a probabilistic estimation problem 8 . 
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Figure 1. Trilateration with and without noise in distance Rj or equally arrival time t\ and/or V g . 

Several approaches have been proposed in literature for damage/impact localization in plate-like structures. Kundu et 
al. 4 proposed a method based on optimizing an objective function to find the impact location in isotropic and anisotropic 
plate. However, the proposed objective function in that reference had the inherent problem of multiple singularities 
which was overcome in Kundu et al. 12 by modifying the objective function. The optimization technique was further 
improved by Ffajzargarbashi et al. 13 . In references 14 ' 15 the nonlinear least square optimization adopting the Gauss- 
Newton method was proposed to determine the location, time lag, and velocity of "synthetic" AE signals. Ciampa and 
Meo 7 used a Newton's iterative method to calculate the coordinates of the impact location and the wave velocity. A 
Genetic Algorithm (GA) for optimization of an objective function based on the modified triangulation methodology was 
proposed in reference 16 . We have previously proposed a reliable probabilistic framework in source localization based on 
Extended Kalman Filter (EKF) estimation 8 . Although the approach improved the accuracy and ability of flexible sensor 
fusion for location estimation by using EKF, it just considered the systematic error in TOF measurement due to 
Heisenberg uncertainty in wavelet method utilized for TOF measurement. And we did not take into account the 
environmental noise effect. Furthermore, implementation of Continues Wavelet Transformation (CWT) in the signal 
processing stage of the algorithm makes it less practical for real time SHM applications, since continues wavelet is in 
general computationally expensive. Here we extend our work to improve the efficiency of the algorithm and consider 
the environmental noise as a source of TOF measurement uncertainty. The paper begins by briefly introducing the TOF 
based algorithm for the AE source location. Section 3 develops a statistical method based on modeled Lamb waves that 
are artificially corrupted by noise to correlate noise level to TOF uncertainty. This is followed by section 4 which 
briefly introduces the EKF theory and its derivation for AE source localization. Section 5 reports the experimental set- 
up followed by the results in Section 6. Finally, the conclusions are given in section 7. 



Proc. of SPIE Vol. 8348 834804-2 



Downloaded from SPIE Digital Library on 05 Jun 201 2 to 1 28.205.27.221 . Terms of Use: http://spiedl.org/terms 



SOURCE LOCATION 



y ,, 



>At I2 .,-'' 



P sensor3 




At,, 


1 — "^--° 


._.— 


sensor4 


---"" 


(x 4 ,y 4 ) 



,^'- '-'-"_"_ _CV7lQsensor5 

sensorlO"" (x,.y/~--^ . ( (^s^s) 

( x i»yJ ~~-0„''sensori 

"° (*,.y,) 



Figure 2. Acoustic emission source and sensors distribution. 

Assuming an arbitrary Cartesian coordinate system, the AE source is at unknown coordinates (x s , y s ) in the plane of the 
plate and the sensors are located at (x„ y t ), as shown in figure 2. It is well known that given the wave velocity of a 
particular Lamb wave mode, at least three sensors are necessary to locate an AE source in plate-like structures. If t m is 
the travel time or TOF to reach the first sensor (master sensor), and At mj are the time difference between master sensor 
and the /th sensor, the following equations can be obtained: 



O m - x s ) 2 + (y m - y s ) 2 - (t m V g ) 2 = 



(1) 



(*! - x s ) 2 + (y t - y s ) 2 - [(t m + At mi )V g f = ( 2 ) 

where V g is the group velocity of the particular Lamb wave mode, which is a function of the product of the frequency/ 
and the plate thickness d: 



V g = F(fd) 



(3) 



"g - r KJ U J 

Knowing the properties of the plate, mass density p, thickness d, Young's modulus E, and Poisson ratio v, V g can be 
calculated by the Rayleigh-Lamb equations 11 . Combining equations (1) and (2) the following equation can be derived: 



At n 



y (xj - x s y + (jj - y s y - yo m - x s y + (y m - y s y 



(4) 



In general, the difference between the time of flight of master sensor to different sensors (At„i) are used to solve this set 
of nonlinear equations with unknowns X=(x S! y s , Vg), However, as mentioned in the previous section, several sources of 
error in the time of flight measurements (i.e. dispersion, noise, temperature, etc.), can affect the accuracy of this 
solution. In this work, to take into account the uncertainties in arrival time t\ or TOF and the wave velocity V g „ the 
arrival time t x and the unknowns x„ y s , V s are treated as mutually independent Gaussian random variables. Based on this 
assumption, the probability density function of the time difference At mj can be also defined as a Gaussian random 
variable 18 with mean and variance defined as: 



At n 



ti 



J At„ 



a t 2 + a t 2 



(5) 



The determination of time differences At mi and the related uncertainties in a noisy environment shall be discussed in 
section 3. In this probabilistic framework, an Extended Kalman Filter (EKF) approach is used for the optimal estimation 
of the state vector X = (x s ,y s , 14). 
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3. STATISTICAL CHARACTERIZATION OF TOF MEASURMENT 

The accuracy of TOF measurement strongly depends on the SNR. In fact, the higher the SNR the lower is the 
uncertainty in TOF measurement. This section carries out a statistical characterization of TOF errors for a small 
ambient noise level. At this aims, first flexural Lamb waves are modeled analytically in order to predict the acoustic 
field in an isotropic plate-like structure generated by an assumed AE source. Then, the simulated signals are corrupted 
artificially with the noise level using a Monte Carlo simulation in which the statistical distribution of the TOF error is 
calculated. 

3.1 AO Lamb wave model as a result of continues point excitation 

In plate-like structures two fundamental Lamb wave modes can propagate: Symmetric (S n ) and Antisymmetric (A n ). 
Displacements of the symmetric modes occur in the direction of wave propagation (i.e. extensional modes), whilst 
antisymmetric modes have displacements transverse to wave propagation direction (i.e. flexural modes). In this paper, it 
is assumed that only the lowest order flexural mode (AO) is the dominant mode. This is a commonly used assumption, 
considering that pencil lead breaks on the surface of thin plate, which have been used during the experimental tests to 
validate the proposed algorithm, predominantly generate the AO Lamb wave mode 19 . Lamb waves have been modeled 
by the Ditri's technique with normal point excitation on the surface of plate 20 . In particular, to model the total surface 
displacements u z associated with an input of finite duration flj), a Fourier Transform is used to decompose the input 
signal into its different frequency components, J{ci>).The input, for pencil lead break simulations, can be assumed as a 
one cycle toneburst with hanning window and center frequency of 1 MHz 21 . The out of plane displacement as a function 
of time and distance from the source is given by 20 : 

u z( r < t) — I /( w ) / E m ((ti)Hl(k m r) exp(— icot) dco 

where t is the time, r is the distance from the point source, k m is wave number of the excited Lamb mode m 
corresponding to co, which can be calculated through the phase velocity V^ (k m = co/V^), and i is imaginary unit; 
E m (a>) is the excitability function for circular crested waves which depend on the angular frequency co and excited mode 
m. The spatial variation in crested circular Lamb wave model is expressed by the Hankel function of first kind Hg. 
Using the method of Viktorov 22 the excitability can be computed. The excitability of a mode is totally governed by the 
amount of out of plane displacement of the mode due to normal point excitation. To model the Lamb wave due to 
preassumed excitation the AO mode is considered as the dominant mode and the effect of all other Lamb modes are 
neglected. Although the assumption of neglecting the effect of other modes is not strictly true, it is almost reasonable 
since most of excitability for normal point excitation and normal displacement field is strongly related to AO 23 . The 
phase velocity V ph for different modes and excitability E(co) just for fundamental antisymmetric AO mode, as a function 
of frequency-thickness are shown in figures 3 and 4 respectively. 
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Figure 3. Phase velocity curves for circular crested Lamb Figure 4. Excitability curves for circular crested Lamb waves 

waves in aluminum plate. in aluminum plate. 

The equation (6) enables the surface displacement to be plotted as function of both distance and time, when the plate is 
excited by a finite duration input at a point. It is assumed that at a sensor position, r, the variation of out of plane 
displacement u z with time is directly proportional to output voltage of the sensor. Figure 5 shows some of the modeled 
signals at different distances to the source. The signals are modeled at range of 5 cm to 50 cm and the duration and 
sampling frequency of all signals is fixed to 0.2 msec and 3 MHz respectively. All the signals are normalized to the 
maximum response at 5 cm. 
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Figure 5. Signal simulation for one cycle sin input in one point with center frequency 1 MHz. 



3.2 



TOF measurement 



Several methods have been proposed in literature for the TOF measurement of acoustic signals, including threshold 
crossing, coloration methods, wavelet transform, and curve fitting approaches" . In this work the TOF has been 
estimated from both experimental and numerical signal using the maximum envelope method 16 . The envelope can be 
calculated using the definition of analytic signal x A (t) of a real signal x(t): 



x A (t) = x(t) + ix H (t) 



(7) 
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Where i is the imaginary unit and x H (t) is the Hilbert transformation of x(t)defined as: 

1 f +co 1 

n]^ t-T 

The magnitude of the analytical signal which is identical to the magnitude of the real signal is called envelope: 



(8) 



A(t) = Vx(t) 2 +x H (t) 2 ( 9 ) 

The time instant in which the envelope of the signal reaches its maximum (A max ), is considered as the TOF. Envelopes 
of some of the modeled signals are depicted in figure 6; also the nominal TOF of the simulated signals for distance of 5, 
10, and 15 cm are shown in figure 6. 
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Figure 6. Nominal TOF for different distances. 
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Figure 7. TOF error definition. 



3.3 TOF error of noisy data 



Investigations on error in TOF measurement of signals due to random noise have been reported for threshold crossing , 
cross correlation 24 ' 26 and some other TOF measurement methods. In this paper the noise modeled as an additive zero 
mean Gaussian random noise with 0.01 standard deviation approximated in the laboratory by means of acquiring several 
signals in 5 cm distance with several pencil lead breaks as AE sources and normalizing the standard deviation of noise 
to the average of maximum signal envelopes. It should be noted that the level of the background noise is measured 
using Root Mean Square (RMS) value of a predetermined portion of signal before the first arrival. In this section the 
TOF error is defined as the deference between nominal TOF captured from modeled signals and TOF extracted from the 
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modeled signals artificially contaminated with the noise. For the sake of clarity figure 7 shows a modeled signal at 50 
cm from the source with and without noise and the corresponding TOF. This figure simply shows a one simulation of 
Monte Carlo for a modeled signal at 50 cm and noise level of 0.01. 
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Figure 8. Monte Carlo analysis. 

To capture the distribution of TOF error and the corresponding statistics, 3000 Monte Carlo simulations were generated 
using the noise level and modeled signals at different distances. From the histograms at the top of the figure 8, it can be 
inferred that as SNR decreases the standard deviation of TOF error (a t ) increases. The results agree with the 
expectation of increasing the TOF uncertainty with respect to distance (decreasing SNR). Figure 9 shows the linear 
relationship between standard deviation of TOF error versus distance to the AE source location. 
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Figure 9. Standard deviation of TOF error versos distance for 0.01 Gaussian random noise and 3000 Monte Carlo analysis. 



4. EXTENDED KALMAN FILTER (EKF) 

In general, a Kalman Filter (KF) is a recursive data processing algorithm that estimates the state of a noisy dynamic 
system 27 . When we talk about the state of a system we mean a vector X consisting of n variables describing some 
interesting properties of the system. To estimate the state a KF processes all available measurements, both accurate and 
inaccurate measurements. The KF has been widely used in a wide range of applications due to its simplicity, optimality, 
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robustness, sensor fusion ability and its capability of filtering out the uncertainties due to the mathematical model and 
measurements ' ' ' . In general the KF is a two steps process: 1) state prediction according to a mathematical model 
and 2) correction of the state according to the measurements collected by sensors. More specifically, in the prediction 
step the KF estimates the state of the system at a given time instant and then obtains a feedback control in the correction 
step, by incorporating a new measurement result into a priori estimate, in order to gain an improved a posteriori 
estimate. Although the underlying approach is very promising there is a critical limitation. In fact the KF assumes that 
system and measurements are adequately modeled by a linear dynamic system. The Extended Kalman Filter (EKF) is 
similar to the KF but it can be used in non-linear systems while making use of first-order Taylor series expansion 
around the estimated mean 27 . For this reason the EKF has been used to handle the set of nonlinear equations (4). 

In this work, the state of the system X consists of three parameters, the AE source location (x s , y s ) and the wave velocity 
V g . Given a sparse array of N sensors, the time differences At mi between the rth sensor and the master sensor represent 
the measurement data vector Z = [At m j] ( - w _ 1 - )xl , N > 4. The state of the system X is related to the measurement vector 
Z through the nonlinear measurement function h defined in equation (4). The key idea underlying the proposed 
approach is that, the AE source location can be considered as unmovable (static) or slowly fluctuating point. Dehghan 
Niri and Salamone 8 here simplified the EKF algorithm for AE source localization as follows: 

z k = Kx^) (io) 

Z fc is predicted measurement vector in kth step according to state estimation X k _ 1 in (k-l)th step using nonlinear 
measurement function h(.) defined in equation 4. 

X k =X k _ 1 + K k {Z-2 k ) (11) 

Where X k is estimated states in kth step derived by means of correcting estimated states at (k-l)th step according to 
difference between actual measurement Z and predicted measurement vector Z k and Kalman gain K k . 



P k = (I- K k H k )P 



fc-1 



Where P k is error covariance matrix in kth step. 



H k = 



dh(X) 



dX 



x=x k . 



The Jacobian matrix H k contains the partial derivatives of the measurement function h with respect to the state X. 

K k — P k-\H kifl kP k-lH k + R) 



(12) 
(13) 

(14) 



/f fc (3x(N-l) matrix) is the Kalman gain used in correction step. 

The objective is to estimate the state vector X = [x s , y s , V g ] given the measurement vector Z and covariance matrix R. 
It should be noted that the measurement vector Z in equation (11) does not have the subscript k because its (N-l) 
components (difference in arrival times) can be considered constant, once the AE signals have been acquired. The 
uncertainty of measurement vector Z is modeled as a zero mean white Gaussian noise with covariance matrix R with 
diagonal terms defined in equation (5). 



R = 



°At ml 







<*lt 



m(n-l) 



(15) 
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The an has been defined according to distance and fixed level of noise (0.01) in section 3.3. Because the location of the 
source is certainly unknown parameter the standard deviation can be selected from the average distance of sensors [28.9 
cm] and then after estimating the location in each step, the R can be modified iteratively. But as a result of low level of 
noise and relatively short distance between sensors the assumption of fixed value for standard deviation of TOF error 
corresponding to average distance is not unreasonable. The standard deviation of TOF error for 28.9 cm and 0.01 noise 
level according to figure 9 is 1.7 x 10~ 7 sec . So from equation (5) <r| t . = 5.78 x 10~ 14 sec 2 . 

The EKF will iteratively correct a priori knowledge of the state vector Xand covariance matrix P, with respect to the Z 
and R. The iteration begins with the initialization of the state vector estimate X and its covariance matrix P : 



Po = 



SO 
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(16) 
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4.1 



Initiation of Location 



The initial location in x and y direction are modeled as an uncorrelated joint Gaussian probability density function. The 
upper and lower bound for each coordinate (for example maximum and minimum coordinates of sensors in x and y 
direction) are computed to calculate the mean and variance of initial location. Then with uniform assumption mean and 
variance can be expressed as: 



_ ( x sU + x sl) 

*s - 2 

_ (y su + y s i) 



(17) 



x^T^su-^sl) 2 



(18) 



Having these means and variances the uncorrelated joint Gaussian probability density function of initial location is 
approximated. Two examples of this joint Gaussian probability density function of initial location are illustrated in 
figure 1 and their parameters are listed in Table 1 . It should be mentioned that by putting the reference of the Cartesian 
coordinate system on sensor 1 the lower bound for x and y coordinate are zero. 



Table 1. Two initial location properties 



Initial Gaussian 
location 


1 


2 


x sU [cm] 


12.7 


25.4 


J s u[cm] 


12.7 


25.4 


a x 2 [cm 2 ] 


13.44 


53.76 


a 2 [cm 2 ] 

y s 


13.44 


53.76 


£s [cm] 


6.35 


12.7 


y so t cm ] 


6.35 


12.7 
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It is noteworthy that this approach is robust against the accuracy of initial guess and its sensitivity will be decreased by 
using the information gathered from more sensors 8 . 



First and second initial Gaussian source location 




Sensor 1 



Sensor 3 



Figure 10. Two different joint Gaussian initial location guesses. 



4.2 wave velocity initiation 



As a result of dispersive behavior of Lamb wave, the uncertainty in frequency corresponding to measured TOF, directly 
causes the uncertainty in group wave velocity assumption even with knowledge of dominant Lamb mode (here A0). In 
this research, we propose to assume that the maximum of Fourier transformation of each acquired signal corresponds to 
its TOF measured by means of maximum envelope approach. It should be insisted that this assumption is not strictly 
accurate since the Lamb wave signals can be totally nonstationary. However, with probabilistic consideration of V g and 
considering it as one of the estimation parameters in this method, its true value will be estimated. If F ; represents the 
maximum of Fourier transformation of ith recorded signal its velocity V gi can be defined form group A0 Lamb wave 
dispersion equation (4) as: 



V gi = F(F t d) 
The mean of V g ; or initial group wave velocity V g is calculated as: 



(19) 



% a 



V 3 



N 



(20) 



And its variance er v 2 easily can be calculated as: 



N 



(21) 
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Examples of this Gaussian probability of V g are illustrated in figures 13 d). 

Signal Processh 




Figure 11. Proposed approach flowchart. 



Overall, the proposed approach consists of two stages as shown in the flowchart in figure 1 1 ; the first stage provides 
information about time of arrivals, frequency and their uncertainties; the second stage uses a priori information about 
the states and iteratively estimate the state vector according to the measured information obtained from the first stage. 

5. EXPERIMENTAL SETUP 

Experiments were carried out using an aluminum plate with dimensions 910mmx910mm><3. 175mm to validate the 
proposed algorithm. The plate was instrumented with an array of four piezoelectric transducers in a square 
configuration, with 254 mm dimension. The experimental set-up is shown in figure 12. For the data acquisition, a 
LeCroy Wave runner oscilloscope with sampling frequency of 3 MHz was used. Each sensor was connected to a 
preamplifier. During testing, AE sources were generated by pencil lead breaks at 4 systematic grid locations. Post- 
processing of the received signals was performed with a PC running a Matlab software code implemented by the 
authors. The number of iteration in EKF stage is fixed to 15. 
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Figure 12. Experimental set up. 
6. RESULTS AND DISCUSSION 



Four AE source were simulated using pencil lead break and the proposed procedure was applied. Figure 1 3 shows the 
signal processing and velocity initiation part of the procedure for point #1 illustrated in figure 14. The results prove the 
validity of the proposed method in fast and accurate location estimation of AE sources. Table 2 listed the actual and 
average estimated position of each source using first and second initial location guesses shown in figure 10. Result of 
2.3 mm error in average in x and y coordinates demonstrates the high accuracy gained with this versatile method. The 
number of iteration to reach these results is 7. Figure 14 depicts actual and estimated source locations. Computational 
time on regular personal computer for each source location using the proposed algorithm took less than 0.03 sec, 
showing this algorithm is appropriate for real-time AE applications using compiled code. 

The results show a slight improvement in accuracy compared to the results achieved in the previous work 8 . Although 
the signal processing stage used to extract time-frequency information in this paper is more sensitive to random noise 
than CWT utilized in the previous investigation 8 , it is significantly more efficient than the CWT. 

Table 2. Actual and estimated source positions 
and estimated wave velocity, (with the first initial 
location) 



Point no: 


1 


2 


3 


4 


x Sa (mm) 


50.8 


101.6 


50.8 


101.6 


x s (mm) 


51.73 


103.4 


50.9 


105.2 


J Sa (mm) 


50.8 


50.8 


101.6 


101.6 


y s (mm) 


51.99 


58.47 


102.0 


104.3 


V g (m/sec) 


3149 


2998 


2891 


2728 


The origin of Cartesian coordinates is at sensor 1 



Using two different initial location guesses which are almost far from each point demonstrates that the fast convergence 
of EKF is highly insensitive to initial guesses unlike other optimization methods like simplex, and Newton. 
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Figure 13. a-d, POINT 1 a) time history of recorded waveforms for four sensors, b) Envelope of waveforms using Hilbert 
transformation, c) FFT of recorded signals d) Initiation of group wave velocity of A0 Lamb mode. 
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Figure 14. Impact location estimation results for four points. 



7. CONCLUSION 



In this paper a probabilistic framework for locating acoustic emission sources is proposed based on the EKF theory. The 
proposed approach, in which TOF and wave velocity are considered as Gaussian random variables, consists of two main 
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stages. In the first stage TOF measurement of Lamb waves are carried out by maximum envelope approach and 
providing the probabilistic information of group wave velocity of AO Lamb mode based on FFT of signals and AO 
group velocity depression curve. Random portion of TOF uncertainty as a result of fixed level of noise are considered 
and modeled numerically by means of modeling signals as a result of predefined normal point excitation and using 
excitability function of AO mode for normal displacement field. The second stage EKF uses the information extracted 
from the first stage and iteratively estimates AE location and wave velocity. In addition to advantages that reported in 
reference 8 like the capability of: (1) taking into account uncertainties in TOF measurements and wave velocity and (2) 
efficiently fuse multi-sensor data to perform an efficient AE source localization, by using more efficient signal 
processing stage in this paper the algorithm becomes more appealing for on-line structural health monitoring. The 
performance of the proposed algorithm has been validated through pencil-lead breaks performed on an aluminum plate 
at systematic grid locations. The plate was instrumented with an array of four piezoelectric transducers. 
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